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We review the recent developments in the field of nuclear double beta decay, which is presently 
an important topic in both nuclear and particle physics. The mechanism of lepton number violation 
within the neutrinoless double beta decay (0^/3/3-decay) is discussed in context of the problem of 
neutrino mixing and the R-parity violating supersymmetric extensions of the Standard model. The 
problem of reliable determination of the nuclear matrix elements governing both two-neutrino and 
neutrinoless modes of the double beta decay is addressed. The validity of different approximation 
schemes in the considered nuclear structure studies is analyzed and the role of the Pauli exclusion 
principle for a correct treatment of nuclear matrix elements is emphasized. The constraints on 
different lepton number violating parameters like effective electron neutrino mass, effective right- 
handed weak interaction parameters, effective Majoron coupling constant and R-parity violating 
SUSY parameters are derived from the best presently available experimental limits on the half life 
of O^/3/3-decay. 
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I. INTRODUCTION 



The double beta {(3(3) decay attracts the attention of both experimentalists and theoreticians for already a long 
period and remains of major importance both for particle and nuclear physics. 

The double beta decay is a second order process of weak interaction and there are few tenths of nuclear systems [0 , 
which offer an opportunity to study it. The (3(3 decay can be observed because the pairing force renders the even-even 
nuclei with even number of protons and neutrons more stable than the odd-odd nuclei with broken pairs. Thus, the 
single beta decay transition from the even-even parent nucleus (A,Z) to the neighboring odd-odd nucleus (A,Z+1) is 
forbidden energetically and the (3(3 decay to the daughter nucleus (A,Z+2) is the only possible decay channel. 

There are different possible modes of the double beta decay, which differ from each other by the light particles 
accompanying the emission of two electrons. We distinguish the double beta decay modes with and without lepton 
number violation. 

The two- neutrino double beta decay (2^/3/?-decay), which involves the emission of two electrons and two antincu- 
trinos, 



is a process fully consistent with the standard model (SM) of electroweak interaction formulated by Glashow M, 
Weinberg H and Salam Q]. This decay mode with obvious lepton number conservation was first considered by 
Mayer in 1935 @. The inverse half-life of 2v(3f3-decay is free of unknown parameters on the particle physics side and 
is expressed as a product of a phase-space factor and the relevant 2z//3/3-decay nuclear matrix element. Thus, the 
measured experimental half lifes of 2z//3/3-decays give us directly the value of the 2z//3/3-decay nuclear matrix elements. 
In this way 2v (3 /3-decay offers a severe test of nuclear structure calculations. 

The 2i//3/3-decay is already well established experimentally for a couple of isotopes. The most favored for the 
experimental study of this rare process is the transition from the ground state + of the initial to the ground state + 
of the final nuclei because of the large energy release. But recently, increased attention is paid also for transitions to 
the 2 + and + excited states of the final nucleus [|6|— p~6| . The half-lifes of the 2i//3/3-decay range from 10 19 up to 10 24 
years [see Table || . Direct counter experiments have observed 2^/3/3-decays for ground state to ground state transitions 
in 48 Ca, 76 Ge, 82 Se, 100 Mo, 116 Cd and 150 Nd. A positive evidence for a 2i//3/3-decay transition to the excited 
state of the final nucleus was observed for 100 Mo []L0| . The geochemical experiments, which observe the double beta 
decay through daughter isotope excesses in natural samples has proved the existence of double beta decay in 82 Se, 
96 Zr, 128 Te and 130 Te. The radiochemical experiment, which observes the accumulation of daughter isotopes under 
laboratory condition, has observed double beta decay in 238 U . The positive signals of geochemical and radiochemical 
experiments are identified with the lepton number conserving 2^/3/3-decay mode. 

The neutrinoless mode of the double beta decay (Oi//3/3-decay) (proposed by Fury in 1939 [Q), which involves the 
emission of two electrons and no neutrinos 



is expected to occur if lepton number conservation is not an exact symmetry of nature and thus is forbidden in the SM 
(Standard Model) of electroweak interaction. The O^/3/3-decay takes place only if the neutrino is a Majorana particle 
(i.e. identical to its own antiparticle) with non-zero mass. 

The study of the Qv(3(3-decsy is stimulated by the development of grand unified theories (GUT's) representing 
generalization the of SU (2) l (g> U(l) SM. In spite of the facts that the SM represents the simplest and most economical 
theory, which describes jointly weak and electromagnetic interactions and in spite of the fact that it has been very 
successful, wherever it has been tested, the SM can not answer many of the fundamental questions, e.g.: Is the neutrino 
really massless? If it has mass, why is this mass much smaller than that of corresponding charged leptons? What 
kind of particles are neutrinos, Dirac or Majorana? Does neutrino mixing take place? Therefore, the SM cannot be 
considered as the ultimative theory of nature. The non-zero neutrino masses and neutrino mixing appear naturally in 
many different variants of GUT's like the simplest SO(10) left-right symmetric model |45| , minimal supersymmetric 
standard model (MSSM) and their extensions. The exception is essentially the minimal SU(5) model which 
however has been practically ruled out due unsuccessful searches for nucleon instability. The expectations arising 
from GUT's are that the conservation laws of the SM may be violated to some small degree. The GUT's offer a 
variety of mechanisms which allow the Oi//3/3-decay, from which the well-known possibility is via the exchange of a 
Majorana neutrino between the two decaying neutrons pMThEg]. 



(A, Z) -> (A, Z + 2) + 2e" + 2v t 



(1) 



{A, Z) -> (A,Z + 2) +2e~, 



(2) 
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If the global symmetry associated with lepton number conservation is broken spontaneously, the models imply the 
existence of a physical Nambu-Goldstone boson, called Majoron |5(i|-|53|], which couples to neutrinos. The Majoron 
might occur in the Majoron mode of the neutrinoless double beta decay f34] , [55| ]: 

(A,Z) -» (A,Z + 2) + 2e- + cf>. (3) 

There are also other possible mechanisms of 0i^/3/?-decay induced by lepton-number violating quark-lepton interactions 
of R-parity non-conserving extensions of the SM J5(],|57| . A complete analysis of this mechanism within the MSSM 
for the case the initial d-quarks are put inside the two initial neutrons (two-nucleon SUSY mode) was carried in [^8| . 
Recently, it has been found that a new contribution of the R-parity violating (lfi p ) supersymmetry (SUSY) to the 
OvftP via pion exchange between the two decaying neutrons dominates over the two-nucleon % SUSY mode |59|. The 
R-parity conserving SUSY mechanisms of Oi//3/3-decay have been proposed and investigated in Ref. p0|] . The GUT's 
predict also new type of gauge bosons called leptoquarks, which can transform quarks into leptons or vice versa ||6l| . 
A new mechanism for 0^/3/9-decay based on leptoquark exchange has been discussed in |]62| . 

The Of/3/3-decay has not been seen experimentally till now. We note that it is easy to distinguish between the 
three decay modes in Eqs. (0)-([|) in the experiment by measuring the sum of electron energies JjJ. A signal from 
the 0i//3/?-decay is expected to be a peak at the end of the electron-electron coincidence spectrums as a function of 
the sum of the energies of the two electrons as they carry the full available kinetic energy for this process. Both, the 
2z;/3/3-decay and the O^/3/30-decay modes yield a continuous electron spectrum, which differ by the position of the 
maximum as different numbers of light particles are present in the final state. The observation of O^/3/3-decay would 
signal physics beyond the SM. It is worthwhile to notice that the 0^/3/?-decay is one of few eminent non-accelerator 
experiments which may probe grand unification scales far beyond present and future accelerator energies. For these 
reasons, considerable experimental effort is being devoted to the study of this decay mode. Presently there are about 
40 experiments searching for the Oi//3/3-decay. The best available experimental limits on the half-life of this decay 
mode for different nuclei are presented in Table [|. The most stringent experimental lower bound for a half-life has 
been measured by the Heidelberg-Moscow collaboration for the Oi//3/3-decay of 76 Ge pl[ : 

T 1 ° / f^- exp (0+ -> 0+) > 1.1 x 10 25 years (90% C.L.) (4) 

The experimental half-life limits allow to extract limits on different lepton number violating parameters, e.g. effective 
neutrino mass, parameters of right-handed currents and parameters of supersymmetric models. 

However, in order to correctly interpret the results of O^/3/3-decay experiments, i.e. to obtain quantitative answers 
for the lepton number violating parameters, the mechanism of nuclear transitions has to be understood. This means 
one has to evaluate the corresponding nuclear matrix elements with high reliability. 

The calculation of the nuclear many-body Green function governing the double beta decay transitions continues to 
be challenging and attracts the specialists of different nuclear models. As each mode of double beta decay requires 
the construction of the same many-body nuclear structure wave functions, the usual strategy has been first to try 
to reproduce the observed 2z//3/3-decay half-lifes. A comparison between theory and experiment for the 2v[3(3 - decay 
provides a measure of confidence in the calculated nuclear wave functions employed for extracting the unknown 
parameters from 0^/3/?-decay life measurements. 

A variety of nuclear techniques have been used in attempts to calculate the nuclear 2z//3/3-decay matrix elements, 
which require the summation over a complete set of intermediate nuclear states. It has been found and it is widely 
accepted that the 2i//3/3-decay matrix elements are strongly suppressed. The nuclear systems which can undergo 
double beta decay are medium and heavy open-shell nuclei with a complicated nuclear structure. The only exception 
is the A=48 system. Within the shell model, which describes well the low-lying states in the initial and final nuclei, 
it is clearly impossible to construct all the needed states of the intermediate nucleus. Therefore, the proton-neutron 
Quasiparticle Random Phase Approximation (pn-QRPA), which is one of the approximation to the many-body prob- 
lem, has developed into one of the most popular methods for calculating nuclear wave functions involved in the j3f3 



decay |63-65|. The QRPA plays a prominent role in the analysis, which are unaccessible to shell model calculations. 
The QRPA has been found successful in revealing the suppression mechanism for the 2i//3/3-decay |63|-|65|] ■ However, 
the predictive power of the QRPA is questionable because of the extreme sensitivity of the calculated 2i//3/3-decay 
matrix elements in the physically acceptable region on the particle-particle strength of the nuclear Hamiltonian. 

This quenching behavior of the 2i//3/3-decay matrix elements is a puzzle and has attracted the attention of many 
theoreticians. The main drawback of the QRPA is the overestimation of the ground state correlations leading to a 
collapse of the QRPA ground state. Several attempts have been made in the past to shift the collapse of the QRPA to 
higher values by proposing different alterations of the QRPA including proton-neutron pairing p6| , higher order RPA 
corrections Q and particle number projection M,p9fl. But, none of them succeeded to avoid collapse of the QRPA 
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and to reduce significantly the dependence of 2^/3/3-decay matrix elements on the particle-particle strength. It is 
because all these methods disregard the main source of the ground state instability which is traced to the assumption 
of the quasiboson approximation (QBA) violating the Pauli exclusion principle. Recently, Toivanen and Suhonen have 
proposed a proton-neutron renormalized QRPA (pn-RQRPA) to study the two neutrino [70[] and Schwieger, Simkovic 
and Faessler the 0^/3/3 decay fn| . In these works the Pauli exclusion principle is taken into account more carefully. 
The pn-RQRPA prevents to build too strong ground state correlations and avoids the collapse of the pn-RQRPA 
solution that plaque the pn-QRPA. The 2^/3/3-decay matrix elements calculated within pn-RQRPA and pn-RQRPA 
with proton-neutron pairing (full-RQRPA) |7lJ has been found significantly less sensitive to the increasing strength 
of the particle-particle interaction in comparison with pn-QRPA results. In the meanwhile, some critical studies 
have discussed the shortcomings of RQRPA like the violation of the Ikeda sum rule |7^,|7J| and particle number 
non-conservation in average in respect to the correlated ground states in the double beta decay J74|]75| ]. 

A large amount of theoretical work has been done to calculate nuclear matrix elements and decay rates. However, 
the mechanism which is leading to the suppression of these matrix elements is still not completely understood. The 
practical calculation always involves some approximations, which make it difficult to obtain an unambiguous decay 
rate. The calculation of the 2z//3/3-decay nuclear transition continues to be subject of interest, which stimulates the 
rapid development of different approaches to the many-body problem, e.g. shell model and QRPA. 

The aim of this review is to present some of the most recent developments in the field of the double beta decay. Let 
us note that several other review articles treating the problem of double beta decay already exist. The theory and the 



experimental status of double beta decay have been reviewed e.g. by Primakoff and Rosen (1981) 76 1, Boehm and 
Vogel (1984) 0, Haxton and Stephenson (1984) §, Schepkin (1984) (47), Doi et al (1985) |§, Vergados (1986) gg, 
Avignone and Brodzinski (1988) @, Grotz and Klapdor ||, Tomoda (1991) f§, and Moe and Vogel (1994) |l|7 

II. MASSIVE NEUTRINOS 

The neutrino remains a puzzle. The neutrino is the only elementary particle, which basic properties are not known 
till today. In contrast with the charged fermions the nature and the masses of the neutrinos has not yet been 
established phenomenologically. The neutrino can be like other fermions a Dirac particle, i.e. is different from its 
antiparticle. However, there is an another possibility proposed long ago by Majorana |32|. The neutrino is the only 
fcrmion which can be a Majorana particle, i.e. it is identical to its antiparticle. This distinction is important if 
the mass of neutrino is non-zero. It is worthwhile to notice that there is no principle, which dictates neutrinos to 
be massless particles as it is postulated in the SM. If the neutrinos are particles with definite mass, it is natural to 
suppose that neutrino mixing does take place. This was first considered by Pontecorvo |8^,Q. It means that flavor 
neutrino fields in the weak interaction Lagrangian are a mixture of fields of neutrinos possessing definite masses. 
Unlike of only one possible scheme for quark mixing there exist several totally different schemes of neutrinos mixing 
as the neutrino is electrically neutral. The number of massive neutrinos depends on the considered scheme of neutrino 
mixing and varies from three to six. The problem of neutrino masses and mixing is very important for understanding 
of fundamental issues o f elem entary particle physics and astrophysics and has been treated in numerous review articles, 
see e.g. |||§||,|||r~ 
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The investigation of neutrino properties is a way to discover new physics beyond the SM. However, the problem of 
neutrino masses and mixing is still far of being solved. The finite mass of neutrinos is related to the problem of lepton 
flavor violation. The SM strictly conserves lepton flavor but the GUT extensions of the SM violate lepton flavor at 
some level. The lepton flavor violation has been discussed in several review articles [ ^9|89| -|9l| . 

By using the flavor fields vl : (^l) c , i^r and {vr) c (the indices L and R refer to the left-handed and right-handed 
chirality states, respectively, and the superscript c refers to the operation of charge conjugation) one can construct 
Dirac, Majorana and Dirac-Majorana neutrino mass Lagrangian's: 

c ° = - M % vil + h - c -> ( 5 ) 

w 

£M = " \ E M n viL+h.c, (6) 

IV 

c d+m = _ [ 1 J^y {m m )vi VlL+ l -vf^ (M%) vl (v m ) c + (M D )in v lL ] + h.c. (7) 

w 

Here the index I and I' take the values e, fj,, r. 
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The Dirac neutrino mass term C D (M D is a complex non-diagonal 3x3 matrix) couples "active" left-handed flavor 
fields vil (weak interaction eigenstates) with the right-handed fields vir which do not enter into the interaction 
Lagrangian of the SM. In the case of the Dirac neutrino mass term the three neutrino flavor field are connected with 
the three Dirac neutrino mass eigenstates as follows: 

3 

vil = ^ Uu v iL . (8) 

i=l 

Here U is a unitary matrix and Vi is the field of the Dirac neutrino with mass m%. The mixing of neutrinos generated 
by Dirac mass term is analogous to the mixing of quarks described by the Cabbibo-Kobayashi-Maskawa matrix. 

It is evident that lepton flavor is violated in the theories with a Dirac mass term, if the matrix M D is not diagonal, i.e. 
additive partial lepton numbers of electronic (L e ), muonic (L M ) and tauonic (L T ) type are not conserved. However, the 
Dirac mass term is invariant in respect to the global gauge transformation vl — > e vl, vr — > e vr (A is arbitrary) 
what implies the conservation of the total lepton number L = L e + + L T . Processes like the Oz//3/3-decay are 
forbidden in this scheme. 

The Majorana mass term C M in Eq. (Q) couples neutrino states of given chirality and its charge conjugate. It 
is obvious that the Majorana mass term C M is not invariant under the global gauge transformation, i.e. does not 
conserve the total lepton number L. Consequently, the theories with C M allow total lepton charge non-conserving 
processes like the O^/3/3-decay. 

The Dirac-Majorana mass term C D+M in Eq. (Q) is the most general neutrino mass term that does not conserve 
the total lepton number L. Mff and M^f are complex non-diagonal symmetrical 3x3 matrices. The flavor neutrino 
fields are superpositions of six Majorana fields Vi with definite masses m, |S6[ , 

6 6 
i=l i=l 

The fields Vi satisfy the Majorana condition 

= v? = C vj, (10) 

where C denotes the charge conjugation and £ is a phase factor. 

The advantage of the Dirac-Majorana mass mixing scheme is that it allows to explain the smallness of the neutrino 
mass within the so called see-saw mechanism proposed by Yanagida, Gell-Mann, Ramond and Slansky J92|. If the 
model involves the violation of lepton number at the large mass scale like in GUT's, the see-saw mechanism is the 
most natural way to obtain light neutrino masses. In the simplest case of one flavor generation the £ D + M takes the 
form 

-D+M _ 1 i 77. \r ( m L m D 



c u+M = ~ ( {y L y vs.) (v L { VR y). (ii) 

2 \rriD tur J 

If we assume the left-handed Majorana mass is zero, the right-handed Majorana mass is of the GUT's scale 
(Mqut ~ 10 12 GeV ) and the Dirac mass mo is of the order of the charged lepton and quark masses (mj? << Mqut), 
we obtain after diagonalization one small and one large eigenvalues: 

mi w F , m 2 « M GUT . (12) 
Mqut 

By considering the general case of three flavor generation we end up with three light masses mi and three very heavy 
masses Mj (m, ~ (m l F ) 2 /Mi). 

The information about the neutrino masses and neutrino mixing is obtained from different experiments. The 
laboratory studies of lepton charge conserving weak processes offer model independent neutrino mass limits that 
follow purely from kinematics. The most stringent bound 4.35 eV on the "electron neutrino mass" has been obtained 
from the Troitsk tritium beta decay experiment |)3| by investigation of the high energy part of the j3 spectrum. 
The tritium beta decay experiments suffer from the anomaly given by the negative average value of m 2 v . It means 
that instead of the expected deficit of the events at the end of the spectrum some excess of events is observed. The 
PSI study of the pion decay gives the strongest upper limit 170 KeV on the "muon neutrino mass" H]. Further 
improvement of this limit is restricted by the uncertainty in the 7r _ mass. By studying the end point of the hadronic 
mass spectrum in the decay r~ — > 27r + 3~^ r an upper "tau neutrino mass" limit of 24 MeV was obtained p4|. 
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The experimental study of neutrino oscillations allows to obtain interesting information about some lepton flavor 
violating parameters. There are neutrino mixing angles and differences of masses-squared Am 2 . In the neutrino 
oscillation experiments one searches for a deficit of some kind of neutrinos at some distance from the source of neutrinos 
vi (disappearance experiment) or one is looking for an appearance of neutrinos of a given kind vy (V = e, /i, r) at 
some distance from the source of neutrinos of different kind vi (I ^ V) (appearance experiment). As a result of 
several oscillation experiments limits have been set on relevant mixing and mass differences. Till now, only the LSND 
experiment p5[ found positive signals in favor of neutrino oscillations. The events can be explained by 7^7 <-> 
oscillations with Am 2 ~ leV 2 . 

The second important indication in favor of neutrino masses and mixing comes from the Kamiokande [ p6"[ , 1MB pTj] 
and Soudan p8| ] experiments on the detection of atmospheric neutrinos. The data of the Kamiokande collaboration 
can be explained by «-> v e or <-> v T oscillations with Am 2 ~ 10~ 2 eV 2 . 



The third most impo rtan t indication comes from solar neutrino experiments (Homestake p9|, Kamiokande [100 



Gallex 1 1 1 ] and SAGE [102]). The neutrinos from the sun are detected by observation of the weak interaction induced 



reactions. The observed events still pose a persisting puzzle being significantly smaller than the values predicted by 



the Standard Solar Model SSM [ L 03 1 . These experimental data can be explained in the framework of MSW [104] 
matter effect for Am 2 ~ 10~ 5 eV 2 or by vacuum oscillations in the case of Am 2 ~ 10~ 10 ey 2 . 

It is worthwhile to mention that all the above mentioned existing indications in favor of neutrino mixing cannot be 



described by any scheme with three massive neutrinos [105]. This fact implies that a scheme of mixing with at least 



four massive neutrinos (that include not only v e , v^, v T but at least one sterile neutrino) has to be considered. 

A prominent role among the neutrino mass experiments plays the Oi^/3/3-decay violating the total lepton number by 
two unit s. A non- vanishing O^/3/3-decay rate requires neutrinos to be Majorana particles, irrespective which mechanism 
is used llOfifl . This theorem has been generalized also for the case of any realistic gauge model with the weak scale 
softly broken SUSY Qg07|. 



The Oi//3/3-decay is a second order process in the theory of weak interaction and neutrino mixing enters into the 
matrix element of the process through the propagator 

< VeL{xx)v T eLtR {x 2 ) >= J2 U^ k U^ k 1 —^ < u k {x x )V k (x 2 ) > ^pC, (13) 

k 

where the propagator of virtual neutrino v k (k = 1, 2...) with Majorana mass m k takes the form 

< M*iP k (**) >=7^—J dV-*(»-*'>^±^. (w) 

(2tt) J V -m% 

Eq. ( p"3] ) can be derived as result of Eq. ([To]) and the following relations: 

VeL = £ U^ kL (C/4 = U ek e~^\e- 2 ^ = (15) 

k 

UeR = Y. U ek^R = Ui k e ia " (16) 

k 

In the case the matrices U L and U R are not independent i.e. in the case of Dirac-Majorana mass term, they obey 
the normalization and orthogonality conditions [p6[ 

2n 2n 2n 

E u lk VvkT = 8w, E < U vk 6 =0, E < Vvk? = (17) 

k=l k=l k=l 

From LEP data it follows from the measurement of the total decay width of Z vector boson that the number of flavor 
neutrinos is equal to three p4[ , i.e. n=3 and I, I' — e, fj,, r. 

If we suppose that the light neutrino (m k <§C few MeV) exchange is the dominant mechanism for the O^/3/3-decay, 
the 0^/3/3-decay amplitude is proportional to the following lepton number violating parameters: 
i) If both current neutrinos are left handed one can separate an effective (weighted average) neutrino mass 

light light 

< m > = E (Uekf 6 m k = E U ek mk. (18) 
k k 
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We note that < m > may be suppressed by a destructive interference between the different contributions in the sum 
ofEq. 



8j) if CP is conserved. In this case the mixing matrix satisfies the condition U e k — U* k where Cfe = ±* is 

Then we have 



the CP parity of the Majorana neutrino [1 



light 

<m v > = ^ \Uek\ 2 Cfe m k- 
k 

ii) If both current neutrinos are of opposite chirality the amplitude is proportional to the factor 

light 

£LR = ^ U ek U e\ & 
k 



(19) 



(20) 



and the amplitude does not explicitly depend on the neutrino mass. We note that in Eqs. (|T^) and ( pp| ) the summation 
is only over light neutrinos. 

If we consider only the heavy neutrino (mi- 3> 1 GeV) exchange O^/3/3-decay mechanism and both current neutrinos 
to be left-handed, we can separate from the Oi//3/3-decay amplitude the parameter 



heavy 



< m„ > 



(21) 



Here, the summation is only over heavy neutrinos. 

Most investigation study the 0^/3/?-decay mechanism generated by the effective neutrino mass parameter given in 
Eq. (|l8|). The value of this parameter can be determined in two ways. 

i) One can extract | < m v > | from the best presently available experimental lower limits on the half-life of the 
O^/3/3-decay after calculating the corresponding nuclear matrix elements. ^From the results of the 76 Ge experiment 
pot one finds | < m v > | < 1.1 eV [ 109 1 . We note that a significant progress in search for O^/3/3-decay is expected 
in the future. The Heidelberg-Moscow and NEMO collaborations are planning to reach sensitivity of 0.1-0.3 eV for 
| < m„ > |. 

ii) One can use the constraints imposed by the results of neutrino oscillation experiments on < m v >. Bilenky 
et al |110| have shown that in a general scheme with three light Majorana neutrinos and a hierarchy of neutrino 
masses (see-saw origin) the results of the first reactor long-baseline experiment CHOOZ [111] and Bugey experiment 
[ 1 1 2 1 imply: | < m„ > | < 3 x lO~ 2 eV. Thus the observation of the Oz^/3/3-decay with a half-life corresponding to 
| < m v > | > 10 eV would be a signal of a non-hierarchical neutrino mass spectrum or another mechanism for the 
violation of lepton number, e.g. right-handed currents [113,126|, $ p SUSY p6|-p9,116| or others [|60p^JTT^|. 



III. DOUBLE BETA DECAY 



In this Section the main formulas relevant for the 2j//?/3-decay, Majorana neutrino exchange and IfL v SUSY mech- 
anisms of O^/3/3-decay are presented. We note that the theory of the 2z//3/?-decay and Majorana neutrino exchange 
mechanism of O^/3/3-decay was discussed in details e.g. in an excellent review of Doi et al [Q. The theory of the 
Ijtp SUSY mechanism of Ou/3/3-decay was presented comprehensively in |58|,|39|]. The readers are referred to these 
publications for the details which are not covered in this review. 



A. Two-neutrino mode 



The 2^/3/3-decay is allowed in the SM, so the effects on the lepton number non-conservation can be neglected. In 
the most popular two nucleon mechanism of the 2v[3(3 - decay process the beta decay Hamiltonian acquires the form: 

H"(x) = ^e(xh a (l - J 5 )Mx) 3*( X ) + h.c, (22) 

where e[x) and v(x) are the field operators of the electron and neutrino, respectively. j a (x) is free left-handed hadronic 
current. 

The matrix element of this process takes the following form: 
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<f\SW\i> =L$1(2*^ L^(p uP2 ,k 1 ,k 2 )J^( Pl ,p 2 ,k 1 ,k 2 ) 

-(pi *-> Pi) - {fa <-> k 2 ) + (pi <-> P2)(fa <-> fc 2 ), ( 23 ) 

with 

J^(pi,P2,Ai,*a) = [ e i ^ +k ^e i ^ +k ^ 



out <Pf\T(J fJl (xi)J„(x 2 ))\p i > in dxidx 2 . (24) 

Here, originates from the lepton currents and J^{x) is the weak charged nuclear hadron current in the Heisenberg 
representation |jS6,117|. p\ and p 2 (k\ and k 2 ) are four- momenta of electrons (antineutrinos), pi and pj are four- 
momenta of the initial and final nucleus. 

The matrix element in Eq. (|2^) contains the contributions from two subsequent nuclear beta deca y pro cesses and 
2i//3/3-decay. They can be separated, if we write the T-product of the two hadron currents as follows ]118[ : 

T(J, i (xi)J v (X2)) = J^{xi)Jy(x 2 ) + O(x 2 o - x w )[J v {x 2 ), J M (xi)]. (25) 

The product of two currents in the r.h.s. of Eq. (|25| ) is associated with two subsequent nuclear beta decay processes, 
which are energetically forbidden for most of the 2z//3/3-decay isotopes. Thus the non-equal-time commutator of the 
two hadron currents corresponds to 2i//3/3-decay process. 

If standard approximations are assumed (only s- wa ves states of emitted electrons are considered, lepton energies 
are replaced with their average value ...) we have [118] 

Jfiv(pi,P2, fa, k 2 ) = -iZM^S^udyk 
x2TTd(E f -Et+ pio + ho + P2Q + fc 20 ), k = 1, 2, 3, (26) 

where, 



M 2 G " T =< 0t| tt / [A k (t/2),A k (-t/2)]dt |0+ >, (27) 



with 



A k (t) = e lHt A k (0)e-* Ht = V ±-L [H[H...[H, A k (0)}...}}. (28) 



71! 

n=0 



Here, |0 + >; and |0 + >/ are the wave functions of the initial and final nuclei with their corresponding energies Ei 
and Ef, respectively. A denotes the average energy A = (Ei — Ef)/2. A k (Q) is the Gamow- Teller transition operator 
A k (0)=Ei T+(ffi) kt k=l,2,3. 

After performing the integration over the time, inserting a complete set of intermediate states |1+ > with eigenener- 
gies E n between the two axial currents in Eq. ( p?] ) and assuming that the nuclear states are eigenstates of the nuclear 
Hamiltonian, one ends up with the well-known form the Gamow- Teller transition matrix element in second order 

2v = v <o;|A fc (0)|l+xl+|A fc (0)|O t +> 

GT 4" E n -Ei + A ■ 1 ' 

So, 2^/3/?-decay can be expressed in terms of single beta decay transitions through virtual intermediate states. The 
form of Mgy in Eq. (^9|) is suitable for approaches using states in the intermediate nucleus (Intermediate Nucleus 
Approach = INA) to 2i//3/3-decay process like QRPA, RQRPA and shell model methods, which construct the spectrum 
of the intermediate nucleus by diagonalization. 

Under the introduced assumptions the inverse half-life for 2^/3/3-decay can be written in the form with factorized 
lepton and nuclear parts: 

[T 1 2 ; 2 (0+ -> 0+)]- 1 = G 2v \M^ T \ 2 , (30) 

where G 2v results from integration over the lepton phase space. For the decays of interest G 2v can be found e.g. in 
refs. Q|4§. 
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It is worthwhile to notice that the form of M^ T in Eq. ( p7| ) as a time integral is very useful for an analytical study 
of the different approximation schemes. It is obvious that the 2^/3/3-decay operator should be at least a two-body 
operator changing two neutrons into two protons. If the commutator of the two time dependent axial currents in Eq. 
( p7| ) is calculated without approximation for a nuclear Hamiltonian consisting of one- and two- body operators, one 
obtains an infinite sum of many-body operators. However, if the single particle Hamiltonian approximation scheme 
is employed, the transition operator is zero |j74|,[75). The QBA and renormalized QBA (RQBA) schemes imply for 
the 2^/?/3-decay transition operator to be a constant. It means that the 2i//3/3-decay is a higher order process in the 
boson expansion of the nuclear Hamiltonian and its higher order boson approximations arc important. This important 
phenomenon can not be seen within INA [ |74| , f75| |. 

In view of the smallness of the nuclear matrix element M"q T in Eq. ( p9| ) the electron £>i/2-wave Coulomb corrections 
to the 2z//3/3-decay amplitude for 0+ s — > 0+ s transition have been examined. It was done first within the closure 



approximation in [119| with the conclusion that the effect is small. The more careful evaluation of the corresponding 



nuclear matrix elements involves the summation over all virtual and 1 states. In the study performed by Krmpotic 



et al [12C ] it is argued this contribution to be significant. If the £>3/2-wave states of outgoing electrons are c onsid ered, 



the nuclear transition can take place through virtual 2~ intermediate states. However, it has been shown in [121] that 
such contribution to the 2i//3/3-decay half-life can be disregarded. 

The inverse half-life of the 2^/3/3-decay transition to the + and 2 + excited states of the final nucleus is given as 
follows: 



V 
with 



[T?; 2 (0+ -. J+)]- 1 = G 3 "(J+) \M% T {J+)\\ (31) 



_ <o;|.4 t (o)|U><uwo)|o+> 

where s — 1 for J = and s = 3 for J = 2. These transitions are disfavored by smaller phase space factors G 2l/ (J + ) 



due to a smaller available energy release for these processes [125]. 



B. Neutrinoless mode 



The presently favored left-right symmetric models of Grand Unifications pioneered by Mohapatra, Pati and Sen- 



janovic [122] are based on a very interesting idea. It is assumed that parity violation is a low energy phenomenon 



and that parity conservation is restored above a certain energy scale. Parity could be violated spontaneously in 



the framework of gauge theories built on the electroweak SU(2)r <S> SU(2)l ® U(l) gauge group [ 123 1 . In addition 



the left-right group structure can be derived from GUT groups SO(10) [124] and the small neutrino masses can be 
explained with help of the see-saw mechanism in a natural way. 

The left-right models contain the known vector bosons (81 GeV) mediating the left-handed weak interaction 
and hypothetical vector bosons responsible for a right-handed weak interaction: 

Wt = cos (W£ + sin QW^ , = - sin £W£ + cos £Wr . (33) 

The "left- and right-handed" vector bosons are mixed if the mass eigenstates are not identical with the weak eigen- 
states. Since we have not seen a right-handed weak interaction, the mass of a heavy vector boson must be considerably 
larger than the mass of the light vector boson, which is responsible mainly for the left-handed weak interaction. 
Within the left-right models the weak beta decay Hamiltonian constructed by Doi, Kotani, Nishiura and Takasugi 



[125 1 takes the form 



nP = 71 [jLM j£+ + KjL " J * + + mR ^ L+ + XjR » J ^ + h - c - ( 34 ) 

Here, k = rj = tan£ and A = (M1/M2) 2 are dimensionless coupling constants for different parts of the right handed 
weak interaction. Mi and M2 are respectively the masses of the vector bosons W\ and W2. The left- and right- 
handed leptonic currents are 

f L (x) = e(x) 7 "(l " 75Ki(z), f R {x) = e{x)r{l + 1 5 Wr(x)- (35) 
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The left- and right- handed current neutrinos {v e L and v e R) are given in Eqs. (|T^) and (|l6|). The nuclear currents 
are assumed to be 



A 



(J°+ (x),J°+ R (x)) = X,^(x-r„)(gvTgAC n) TgA<T„ + gvD n ), (36) 

n=l 

where C n and D ra are nucleon recoil terms 0] . The r\ term arises when the chiralities of the quark hadronic current 
match those of the leptonic current, i.e. there are of left-right combination (this is possible due W-boson mixing). 
The A term arises when both the quark hadronic current and leptonic current are right-handed [ see Fig. |lj]. 

Let suppose that Oz//3/3-decay is generated by the weak interaction Hamiltonian in Eq. (]34|) and that only exchange 
of light neutrinos is considered. Then, the matrix element for this process depends on the three effective lepton 
number violating parameters < m v > (defined in Eq. (|l8|)), < A > and < r\ >: 

light light 

< A >= A e LR = \ J2 U ek U," k < v >= V £LR = V J2 U ek U ek (37) 

k k 

We note that the Oz^/3/3-decay amplitude converges to zero in the limit of zero neutrino mass. It is apparent by 
< m u >= in this case. If all neutrinos are massless or are massive but light , < A > and < 77 > vanish due to the 
orthogonality condition in Eq. ( |l7| ) . 

The inverse half-life for O^/3/3-decay expressed in terms of the lepton number violating parameters is given by [48[] 

[Tffa]- 1 = \M% T \ 2 [C mm { — ) 2 + C\\ < A > 2 +C m < V > 2 + 

, A <m v > \ T „ , , I < m v > I 

C m x\ < A> y- ^COS*i + C m ^ <T)> N ^cos* 2 + 

m e m e 

C Xv \ < A> || <7y> |cos(*i-* 2 )], (38) 

where and ^2 are the relative phases between < to„ > and < A > and < rj >, respectively. The coefficients C xy 
are expressed as combination of 8 nuclear matrix elements and 9 kinematical factors: 

C mm = (1 — xf) G i, C mv = (1 — Xf)[x2+Gq3 — Xi-Goi — XpGq5 + XrGog], 

1 2 

C m \ = —(1 — Xf)[X2-G 3 — X1+G04], C\x — [X2-G()2 + gXl-G 4 — gXl+X2-G 0a ] 



1 2 „ 2 

1, J 



Grpj — [X2+G02 + -X1-G04 H — prXi-X2+G 3 + XpG 8 — XpXrGq7 + X 2 rG q] 



with 



Cm, = -2[x2-X2+G 02 - - (X1+X2+ + X2-Xi-)G 3 + gXi+Xi-G 4], (39) 



1 . . 

Xl± = ±3XF' + XGT' - 6% T ', X2± = ±XFuj + XGTu - TTXli- (40) 



The numerical values of the kinematical factors Goi (i = 1,9) can be found e.g. in p8,126|. The nuclear matrix 
elements within closure approximation are of the form 

Xj = {Qf\J^T+T+O I (T i ,T j ,<r i ,<r j )\0+) (41) 

with X[ = M\} v (I = GT) or X/ = X i (I = GT, F, GTui, Fuj, GT , F', T', P, R). The radial and spin-angular 
dependence of the transitions operators Oi is given in Table [fj. We mention that the inaccuracy coming from the 
closure approximation in completing the sum over virtual nuclear states is small (about 15%) |127|, since the average 
neutrino momentum is large compared to the average nuclear excitation energy |79[ ]. The nuclear matrix elements 
derived without closure approximation can be found in [128 - 131 1 . The expressions for the O^/3/3-decay matrix elements 
obtained within the relativistic quark confinement model are given in [ 132 1 . 

We note that a more general expression for the inverse half-life including also the heavy neutrino exchange mech- 
anism can be found in [129]. We shall discuss this type of Oi^/3/3-decay mechanism in context of the Ifip -SUSY 
mechanisms. 
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C. Majoron-Neutrinoless mode 



In many extensions of the standard model |5C , 132 , 134 , 54| appears a physical Goldstone boson called Majoron, 
which is a light or massless boson with a very tiny coupling to neutrinos. The Majoron offers a new possibility for 
looking for a signal of new physics in double beta decay experiments. 

The Majoron emitting mode of double beta decay (Of/3/30-decay) [see Eq. (J3J) ] has been proposed by Georgi, 
Glashow and Nussinov |32| within the model introduced by Gelmini and Roncadelli JjlJ . 

Today, we know several Majoron models based on different motivations leading to a different form of the Majoron- 
neutrino interaction Lagrangian. The classical Majoron Models [50,|l 135 suggest that the interaction of the so called 
ordinary Majoron with neutrinos is of the form: 



C<t>w = ~o ^ v i{a>ijPL + bijP R )vj<j)* + h.c. 



(42) 



where Pr,l = 1/2(1 i 75). If the usual approximation rrii.j <C q 
assumed, one can write for the inverse half-life of the 0z//3/3(/>-decay 



Pf ~ O(100MeV) valid for light neutrinos is 



[Tfff]- 1 = \<g> \ 2 \M°J T \ 2 \1 - Xf\ 2 G b , 



(43) 



with the effective Majoron coupling constant 



< 9 >= 



light 



ij 



(44) 



We see that O^/3/30-decay is ruled by the same nuclear matrix elements as the light Majorana neutrino mass mechanism 
of the 0^/9/?-decay. The numerical value of the kincmatical factors Gb for different nuclei can be found e.g. in E§]. 



Recently, a new class of Majoron models [136-13S] have been considered, which introduce a new types of Majorons 
carrying leptonic charge. The theoretical analysis of the charged Majoron modes of double beta decay have been 
performed in [|55,139|. It has been found that within the new Majoron models unobservable small decay rates for the 
Majoron-emitting Qv (3 (5(f)- decay are expected. 



D. SUSY neutrinoless mode 



Supersymmetry has been introduced to solve the so called "hierarchy problem" of the unified theories ]140[ . It is the 
only known symmetry which can stabilize the elementary Higgs boson mass with respect to otherwise uncontrollable 
radiative corrections. Unlike the symmetries of ordinary particle physics, that relate the particles of the same spin, the 
SUSY is relating bosons to fermions and vice- versa and requires the existence of new supersymmetric particles. To each 
ordinary particle like leptons (neutrino, electron), quarks, gauge bosons (W, Z, photon), Higgs bosons corresponds a 
superpartner with different spin (bosons replacing fermions and vice versa), i.e. slepton (sneutrino, selectron), squark, 
gaugino (wino, zino, photino), higgsino. SUSY is not a good symmetry in the nature. Supersymmetry has to be 
broken otherwise the superpartners would possess the same masses as their ordinary partners. This would change 
drastically the phenomenology of particle physics. 

The simplest SUSY extension of the SM is the minimal supersymmetric standard model (MSSM) . The Lagrangian 
of the MSSM conserves a new quantity called R-parity: 



(45) 



Here, B, L and S are respectively the baryon number, lepton number and spin. The R-parity has been imposed guar- 
antee the baryon and lepton number conservation. Ordinary particles have R-parity even and their supersymmetric 
partners R-parity odd. However, R-parity conservation is not required by gauge invariance or supersymmetry and 
might be broken explicitly or spontaneously at the Planck scale ($ p SUSY) [134]. 

The Ifip SUSY is an other way of the lepton number violation in addition to the Majorana neutrino mass term. It 
has interesting phenomenological consequences. The lepton flavor violation in the context of SUSY models has been 
discussed extensively in the literature [14S- 152] . The collider experiments |141| ], low-energy processes [142], matter 
stability [143-145 and cosmology [146- 148 1 offer important constraints on ^ p SUSY theories. 
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The O^/3/3-decay induced by exchange of superparticles has been studied within the lfi p MSSM [^8 59 116]. The 
relevant ]fi p violating part of the superpotential takes the form: 



= \' ijk Li Qj D k . 



(46) 



Here L and Q are respectively lepton and quark doublet superfields while D denotes down quark singlet superfields. 
A^- fc is a coupling constant and indices k denote generations. 

The nuclear 0!//3/3-decay is triggered by the Qv(3f3 quark transition d + d — > u + u + 2e~ . If this lepton number 
violating process (AL e = 2) is induced by heavy SUSY particles and heavy neutrinos in a virtual intermediate state, 
one can write for the corresponding effective quark-electron Lagrangian [ 153 1 : 



C 



AL C =2 
eff 



2m 



e(l + 75 )e c 



VPS JpsJpS " ~£]T Jt ' Jthv + VN JvA^VAfi 



(47) 



The color singlet hadronic currents are Jps — u a (l + 75)^0, = u a a^ v {\ + -f5)d a , J^ v — u a j fl (l — "fb)d a , where 
a is a color index and a^ v — (i/2)[7 A1 , r y v \. Compared to Ref. |38| the Lagrangian in Eq. ( [47| ) contains properly taken 
into account contribution of the color octet currents. 

The effective lepton-number violating parameters VPS and vt in Eq. 
% MSSM as follows: 



fl) accumulate fundamental parameters of 



with 



VPS = V X e + V x f + Vx + Vg + 7v' S , 
Vt =V X - T l x f+Vg-V' g , 



(48) 
(49) 



ira s \'i U mp 



1)9 



Vx 



, F m iR m~ g 



na 2 X11 



i ji/j ; L — ' to 

b d R i=l 



mp 

Xi 



b d R v L ' L 



4 



'xf 



ira s Ajjj m P 
12 G|m? m 5 

b d R » 



W) + e 2 Li (u) 
4 4 



2 / 



2 



7T«2 



2 G, 



S \ m er / TO V . 



£i?i(rf)eLi(e)+ 



+ eLj(M)effi(d) 



+ £Li{u)e Ll (e) 



(50) 
(51) 
(52) 
(53) 
(54) 



where a 2 — g|/(47r) and a s = g|/(47r) are SU(2) L and ST/(3) C gauge coupling constants, to Si _, «i<? an d m x 

are masses of the u-squark, d-squark, gluino g and of the lightest neutralino Xj respectively. The neutralino is linear 
superposition of the gaugino and higgsino fields: \ — a x B + (3 X W 3 + 5 X H® + J X H 2 . Here W 3 and B are neutral 
SU(2)l and U(l) gauginos while H®, are higgsinos which are the superpartners of the two neutral Higgs boson 
fields H® and H® with a weak hypercharge Y = —1, +1, respectively. The mixing coefficients a x , /3 X , 7x , <5 V can be 



obtained from diagonalization of the 4x4 neutralino mass matrix 154]. Neutralino couplings are defined as [154] 
— —T^{iji)P x + tan9w (Ts^ip) — Q(ip))ct x , tp^p — Q(ip) tan9wo: x . Here Q and T3 are the electric charge and 
the weak isospin of the fields ip = u,d,e. The relevant Feynman diagrams associated with gluino and neutralino 
contributions to the Qv(3(3-decay are drawn in Fig. ^|. 

The effective lepton-number violating parameters vn hi Ec[- ( [47| ) has origin in the exchange of heavy Majorana 
neutrino and takes the form: t]n — m p < m~ l >. Here, m p is the mass of proton and < to" 1 > is the inverse value 
of the effective heavy Majorana neutrino mass defined in Eq. (|2l|). 

For the calculation of the O^/3/3-decay nuclear transitions the quark-lepton interaction in Eq. p7| ) has to be 
reformulated in terms of effective hadron-lepton interaction. The quarks undergoing ]fi p SUSY transition can be 
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incorporated in two nucleons (the well-known two-nucleon mode), one intermediate pion and one nucleon (so called 
one-pion mode) or in two intermediate pions (so called two-pion mode). These three possibilities are drawn in Fig. 
||. We note that intermediate SUSY partners are heavy particles and thus in the two-nucleon mode the two-decaying 
neutrons must come very close to each other. This contribution is therefore suppressed by the short range nucleon- 
nucleon repulsion. Therefore, it is expected that the pion-exchange mechanism dominates over the conventional 
two-nucleon one. The scale on which the pion contribution is enhanced compared to the two-nucleon mode is the 
ratio of the nucleon form factor cut-off parameter compared to the pion mass. In addition, the relative importance 
of the two-nucleon mode, of the one-pion mode and of the two-pion mode depends on the evaluation of the hadronic 
matrix element < 7r + |JiJi|7r _ > (i = P,S,T), which is model dependent. The values of the corresponding structure 
coefficients calculated within the vacuum insertion approximation (VIA) and within the non-relativistic quark model 
(QM) differs from each other by a factor of three |)9| . It is worthwhile to notice that the importance of the pion- 



exchange currents in the CV/3/3-decay was first been pointed out by Pontecorvo [155]. 

The half-life for the neutrinoless double beta decay regarding all the three possibilities of hadronization of the 
quarks can be written in the form |153| 



G 



01 



[z?; a (o+-o+)] 

VrMf + (tips - r] T )M 



2N 
f 



+ i(VT + ivPs)M 7TN + t] N M 
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(55) 
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(56) 



Here Goi is the standard phase space factor introduced in Eq. ( |39| ) and ttla — 850 MeV is the nucleon form factor 
cut-off (for all nucleon form factors the dipole shape with the same cut-off is considered). The nucleon structure 
coefficients a's entering nuclear matrix elements of the two-nucleon mode are given in Table III. The structure 



coefficient of the one-pion a 1 * and two-pion mode a 27r are §9,153]: a l7r = -0.044 and a 2n = 0.2 (VIA), 0.64 (QM) 



The partial nuclear matrix elements of the ]ft, p SUSY mechanism for the Ov (5 '/3-decay appearing in Eq. ([56|) are of the 
following form: 



Mi = (0+| £ t+t+Oj (r< , r, , a t , a, ) \0j 



(57) 



where I = FN, GTN, F', GT , T', GT - in and T - in (i=l,2). The corresponding operators 0/ are given in Table 
After the nuclear matrix elements are evaluated one can extract the constraints on the combination of lepton number 



violating parameters r]ps, T)t and tjn from the non-observation of the Oi//3/3-decay [see Eq. (|55j)]. It is worthwhile to 
notice that it is a common practice to consider separatelly $ p SUSY and heavy neutrino exchange mechanisms for 
the 0^/3/9-decay. Then, the experimental lower bound T^J£ f° r t ne half-life provides the following constraint on t/n 
and r] SUSY parameters: 



l]N 



< Vn P = 



\M 



N 



Ct oiJ 1/2 



(58) 



»7 S! 



3 V 



\M 



vN\ 



1 



G m F 



cxp 
1/2 



(59) 



where we assumed the dominance of the two-pion mode (, which seems to be very probable). Then we have for the 1st 
generation Yukawa coupling constant Am under the phenomenologically viable simplifying assumptions (the lightest 
neutralino is B-ino dominant, rriq > m§/2) [153]: 
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X; il <12.5^( i5 ^) 2 (i) 1/2 . (61) 

We note that the running QCD coupling constant a s (Q) was taken at the scale Q — lGeV with the normalization on 
the world average value a s (M z ) = 0.120 Q. 

There is also a possibility to get useful constraints on other parameters of the SUSY models by (e.g. B-L violating 
sneutrino mass) considering the R-parity conserving SUSY mechanism for O^/3/3-decay. A class of such mechanisms 
have been considered in [ pOfl . However, the constraint on the sneutrino mass spectrum from the double beta decay 
experiment has been found less stringent then the constraint from the best laboratory limits on the neutrino masses. 



IV. NUCLEAR MODELS 



A. Shell model 



The shell model is giving a very satisfactory description of the observed properties of the low-lying states of nuclei. 
The microscopic shell model wave functions posses a set of quantum numbers reflecting the symmetries of the strong 
interaction. This is an advantage in respect to the QRPA (Quasiparticle Random Phase Approximation), where the 
proton and neutron numbers are only conserved in average. However, for the medium and heavy nuclei the shell model 
space increases so drastically, that reliable calculations are impossible beyond the pf-shell. But even relative large 
spaces in pf-shell nuclei can not reliably describe the states in the giant Gamow- Teller resonance region, where most 
of the (3~ strength is concentrated. Thus even in lighter nuclei the shell model studies of the 2^/3/3-decay transitions 
are based on the assumption that the contribution of higher lying (3~ strengths to 2i//3/?-decay matrix element is 
negligible, e.g. due to the small overlap with the (3 + strength distribution and due to the large energy denominator. 

Calculation of the O^/3/3-decay matrix elements requires only the knowledge of the initial and final nuclear wave 
functions as the exchange potential depends weakly on details of nuclear structure. It is worthwhile to notice that 
it is almost impossible to construct even the wave functions of the initial and final nuclei for nuclear systems, which 
are interesting for the (3(3 decay. Only for the double beta decay candidate 48 Ca this has been found possible by 
considering the full (pf) shell [156|. In the heavier nuclei the number of basis states increases astronomically and it is 
difficult to perform a realistic shell- model calculation without very severe truncations. 

In earlier studies the weak coupling scheme has been employed, starting from a product of neutron and proton 
wave functions and employing truncations according the energies of the unpertubated proton and neutron states 
Jl,157 15S | . Nowadays, there are some shell model codes [159-162] which can treat some of the heavier nuclear 
systems undergoing double beta decay in a considerably more realistic way. We remind also that the feasibility of the 
shell model calculations is growing with increasing computer facilities. 

The shell model calculations generally have resulted in larger transition matrix elements than the QRPA. One 
may surmise that the stringent truncation imposed on shell model spaces exclude configurations responsible for the 
redu ction by des truc tive interference. This idea is supported by the shell model studies of the 2zy/3/3-decay of 48 Ca 
[ 160 1 and 92 Mo [161|, in which the dependence of the 2v (3 /3-decay matrix elements on the configuration space used 
have been examined. Another reason can be the omission of the spin-orbit partners. Recently, it has been found that 
there is a destructive interference between "spin- flip" and "non spin- flip" contributions to the 2i//3/?-decay transition 
leading to the suppression of 2i//3/3-decay matrix elements [164,165]. In addition the matrix elements obtained within 
the shell model can be suppressed by the renormalization of the Gamow- Teller operators in order to account for 
about 20% quenching of the beta decay strength. Till now, the origin of this quenching phenomenon is not well 
understood and can be explained by a number of reasons, e.g. by influence of many-particle configurations or by 
the renormalization of axial vector and vector weak interactions constants due to nuclear medium effects leading to 
9a/9v ~ 1. 



B. QRPA and renormalized QRPA 



The pn-QRPA has been found to be a powerful model, considering its simplicity, to describe nuclear matrix elements 
in ordinary beta decay. In addition, the QRPA can handle of great number of intermediate states. Therefore, the 
pn-QRPA has been the most frequently used nuclear structure method to deal with the nuclear structure aspects of 
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the Pf3 decay for open shell systems [see e.g. |63| |65|Jl29| , |l30|Jl66| - |i68| ]. In the double beta decay the initial nucleus 
decays to final nucleus through virtual excitations of all states of the intermediate nucleus. The excited states are 
obtained by solving the well-known QRPA equations, which yield all excited states of a given even-even nucleus. The 
advantage of the QRPA calculation in respect to the shell model approach is that it can be performed in large model 
spaces. The QRPA produces a considerably smaller number of states than the shell model as it is based on a harmonic 
approximation. 

A shortcoming of the pn-QRPA is that it frequently overestimates the ground state correlations, leading to insta- 
bilities and finally to collapse of the QRPA solution. The collapse is caused by generation of too many correlations 
with increasing strength of the attractive proton-neutron interaction. Unfortunately, the physical value of this force 
is usually close to the point in which the QRPA solutions collapse. This behavior has its origin in the quasi-boson 
approximation (QBA) violating the Pauli exclusion principle and causing the QRPA excitation operators behave like 
bosons |7^,^l|. The increasing violation of the Pauli principle is generated by excessive ground state correlations. 
To overcome this difficulty the renormalized QRPA have been proposed [ |70|]7l| l, which take into account the Pauli 
exclusion principle in an approximate way. The price paid for bypassing the collapse in the RQRPA is the violation 
of the Ikeda sum rule |f^ , ^| . 

The pn-QRPA full-QRPA |§, pn-RQRPA (7(} and full-RQRPA @ formalisms are fairly well known. So 

we briefly mention only its main features here. 

The QRPA as well as the RQRPA formalisms consist of two main steps: (i) The Bogoliubov transformation smears 
out the nuclear Fermi surface over a relatively large number of orbitals (ii) The equation of motion in the quasiparticle 
basis determines then the excited states. 

If the proton-proton, neutron-neutron and proton-neutron pairing is considered, the particle {c± am and c Tam , 
t = p,n) and quasiparticle {a~^ ama and a nam ^ fj, = 1,2) creation and annihilation operators for spherical shell model 
states (The label a designates quantum numbers n a , l a ,j a .) are related each to other by the Hartree-Fock-Bogoliubov 
(HFB) transformation: 
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(62) 



where the tilde ~ indicates time reversal (a TQ 



i-iy 



). The occupation amplitudes u and v and the 



single quasiparticle energies E aa are obtained by solving the adequate HFB equation 



considered is finite, the pairing interactions are r enormalized by the strength parameters d pp , u»„„ 



empirical gaps defined by Moeller and Nix 169 1 . The numerical values of the d pp 



Since the model space 
d 

and d nn are close to unity. A d pn 
value higher than unity is the price paid for the spherical symmetry of the model which excludes the treatment of the 
T=0 pairing. The J=0 T=0 pairs can be treated in a BCS or even HFB approach only due to deformation. The T=0 
pairing is effectively taken into account by the renormalization of the T=l J=0 n-p interaction leading to a higher 



value of d pn . For large enough model space the value of d pn is about 1.5, what is reasonable 110' 



In the limit without proton-neutron pairing one has: u 



2p 



J 2p 



transformation in Eq. ( |62| ) is reduced to two BCS transformations, for protons (u lp 



0. In this case the Bogoliubov 
= u„, Ui„ = i> ) and for neutrons 



(u 2 



v n ) separately. 



In the framework of the full-QRPA or full-RQRPA the m th excited states with angular momentum J and projection 
M is created by a phonon-operator Q with the properties 



Q 



JM\"RPAI 



\m,JM) 



and 



Q|0+ 



RPAI 



0. 



(63) 



Here, |0pp is the ground state of the initial or the final nucleus. The phonon-operator Q takes the following form 



Q7L = E X i™( fc < J )4 2 (fc, I, J, M) + Y$(k, I, J)A 12 (k, I, J, M) + 

kl 

E X ™^ k > l > J ) A U k > l > J < M ) + O fc > J) A ^(k, I, J, M). 



(64) 



k<l 



Aj lu (k, I, J, M) is two quasiparticle creation and annihilation operator coupled to angular momentum J with projection 
M namely 
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JM a f a 1 



Al v (k,l,J,M) = n(kfi,W) C 'i 

m k ,mi 

n(fc/x, M = (1 + (-1) J *« + hi W 3/2 - 



(65) 



We note that if proton-neutron pairing is neglected, the phonon operator in Eq. (|64|) decouples on two phonon 
operators, one for charge changing and second for non-charge changing modes of nuclear excitations. 

In the full-RQRPA the commutator of two bifermion operators is replaced with its expectation value in the correlated 
QRPA ground state \0qrpa > (renormalized quasiboson approximation). We have 



(Qqrpa\ J- J, M )> A t'»'( k '> V > J ' M )] Krpa) 

j 1 _ J { < {) QRPA\\ a ll^l\^QR.PA 



jk+jl — J 
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(66) 



with jk = \/2jk + 1- If we replace \0q RPA > in Eq. (|6(]) with the uncorrelated HFB ground state, we obtain the 
quasiboson approximation (i.e. V ^k,vi;J" = 1); which violates the Pauli exclusion principle by neglecting the terms 
coming from the Fermion commutation rules of the quasi-particles. The QBA assumes that pairs of quasiparticles 
obey the commutation relations of bosons. 



The full-RQRPA takes into account the Pauli exclusion principle more carefully. From Eqs. (63) and (|6J) one can 
derive the RQRPA equation 
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(67) 



A, B 



x m J m 



X ,Y 



The matrices A and B are given as follows: 

Kpa\ [A a p(a,b,J,M), [H,Al u (k,l,J,M)]] \0+ PA ) 



a aabf3 



B 



aabft 
kfJ,lp,,J 
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RPA\ 



A a/3 (a,b,J,M), HtA^fal, J,M) 







RPAl i 
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(68) 
(69) 



where H is the quasiparticle nuclear Hamiltonian ]7l| . The coefficients T> tl k,vi\J 1 ' are determined by solving numerically 
the system of equations @,fnj : 



V, 



' , I 2 



4e ^v^i ,2 in"'(M',j /7r ')i 



(70) 



The selfconsistent scheme of the calculation of the forward- (backward-) going free variational amplitude X (Y ), 
the energies of the excited states for every multipolarity J T and the coefficients T> fl k,vl;J'" ^ s a double iterative 
problem which requires the solution of coupled Eqs. (67) and (f70|). 

In order to calculate double beta transitions two fully independent RQRPA calculations are needed, one to describe 
the beta transition from the initial to the intermediate nucleus and another one for the beta transition from the 
intermediate to the final nucleus. For that purpose the one-body transition densities of the charge changing operator 
has to be evaluated. They take the following form: 



< J^rm || [c+ k Cni]j || 0+ >= V2J + 1 ^ m(jik,vl)x 



(71) 



16 



< 0/ II Hk'Zni>]j II J * m f >= V2J+T ^ m(jik\ vl 1 ) x (72) 

H,v=X,2 



with m(fm, vb) = ^^j,^ sly/z ■ We note that the Jt^ v (k,l, J v ) and Y^ v (k,l, J 71 ') amplitudes are calculated by the 
renormalized QRPA equation only for the configurations \xa < vb ( i.e., \i \x — v the orbitals a and b are ordered 
a < b, and if /a = 1, i/ = 2, the orbitals need not to be ordered). For the amplitudes X „ v {k, I, J w ) and F (fc, 7, J 7r ) 
in Eqs. ([n]) and j7^ ) exist the following symmetry relations: 

x; v (k,i,j*) = -(-iy^'- T x^(i,k,r), Y™(k,i : r) = -(-iy^'-'Y^(i,k,r). (73) 

The index i (f) indicates that the quasiparticles and the excited states of the nucleus are defined with respect to the 
initial (final) nuclear ground state |0* > (|0i >). We note that for V^ivj^ — 1 (i-e. there is no renormalization) and 
U>2p — v 2p = Ui n = v\ n = (i.e. there is no proton-neutron pairing), Eqs. ( ffl| ) and (|7^) reduce to the expressions of 
the standard pn-QRPA. 

The two sets of intermediate nuclear states generated from the initial and final ground states are not identical in 
the QBA or RQBA schemes considered. Therefore the overlap factor of these states is introduced in the theory as 
follows [[75): 

< Jm^m, >~ \Q?M^QjM ] ~ E hk'h'Ukp„'U lv V X (74) 

fik<ul, \i' k' <~v' V 

(X: (k, i, TpQi, (k, i, j) r; n ; (k, 1, j)T;,i, (k, i, j)) , 

with 

fiw=4W2+4W2- (75) 

Here, Q™/f and Q^^}' are respectively phonon annihilation and creation operators for the initial and final nuclear 
states. We note that in the previous calculations u was approximated by unity. However, in that case the overlap 
factor is dependent on the phases of the BCS/HFB occupation amplitudes u's and v's, which are in principal arbitrary. 
A negligible difference between the results with the above overlap containing u and with the overlap with u = 1 one 
obtains only if the phases of the BCS/HFB amplitudes are chosen so that u is positive for each level. 

The double beta decay nuclear matrix elements are usually transformed to ones containing two-body matrix elements 
by using the second quantization formalism. One arrives at the expression: 



<0i 2 >= v {-r + ^ +j+j {2j+i)\ Jk 31 J \x 

I Jl' Jk' J \ 

klkl v •* 

< 0/ II [c + pk ,~c n v]j || rm s >< r ms \r mi >< r mi \\ [c+ k c nl \j \\ o+ > 

x <pk, P k l -J\f{r 12 )T+T+0 12 f{r 12 )\nlM'\J> ■ (76) 

with short-range correlation functions 

f(r) = 1 - e- ar \l - br 2 ) with a = l.lfm 2 and b = 0.68fm 2 , (77) 

which takes into account the short range repulsion of the nucleons. In the case of the nuclear matrix element M"q T in 
Eq. (H) 0i2 is replaced by 2a\ ■ a 2 /{^+ + fi™ 7 ) and in the case of neutrinoless double beta decay matrix elements 
in Eqs. and (||) 12 is replaced by (9j(n,r 2 , ai, cr 2 ) [see Table ||. 

In the QRPA and RQRPA calculations it is necessary to introduce the renormalization parameters g pp and g p h 
for particle-particle and particle-hole channels of the G-matrix interaction of the nuclear Hamiltonian H, which in 
principle should be close to unity. Slight deviations from unity can be expected because the full Hilbert space can 
not be taken into account. The g p h value is usually fixed by physical arguments, e.g. by fitting the position of giant 
Gamow- Teller resonance. The g pp is considered as free parameter and the dependence of the matrix element on g pp is 
carried out. The problem regarding g pp is crucial to predictability of the QRPA calculation. It has been shown that 
(3 + transitions (second branch of the double beta decay in reverse) to low- lying states are very sensitive to g pp . The 
introduction of g pp as a phenomenological factor makes it possible to reproduce the observed 2z//?/3-decay half-lifes. 
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C. Operator expansion method 



The Operator expansion method (OEM) is a nuclear structure method for the 2^/3/3-decay, which has the advantage 
to avoid the explicit sum over the intermediate nuclear states. 

Th ere are two different ways to derive the 2z//3/?-decay OEM transition operators. In the first approach (OEM1) 
|{L7l| the expa nsion of the denominators of M^? T in Eq. (29) in Taylor series is used. In the second approach proposed 
(OEM2) [172 1 the OEM is derived from the integral representation of the nuclear matrix element M^ T in Eqs. fl27| ) 
and (|28|). It has been found that the OEM2 offers advantages over OEM1 as there are no problems with convergence 
of the power series expansion of the denominator |74|,[75|], which has been a subject of criticism 173]. In addition, the 
approximations under consideration are better controlled. 



The OEM is based on two main assumptions 171 jlTS ] : i) It is assumed that the kinetic energy operator can be 
ignored in the resulting commutators [ see Eqs. (|27|) and (28|) ] and therefore the nuclear Hamiltonian H is represented 
only by two-body interaction terms, ii) Only two-body terms are retained in evaluating each commutator and higher 
order terms are neglected. These approximations are supported by the following arguments: First, if we consider 
only the one-body part of the nuclear Hamiltonian Hq, the 2z//3/3-decay transition operator [A(t/2), A(-t/2)] in j27| ) 
is an one body operator [ [74|j75| . As the 2^/3/3-decay operator should be at least a two-body operator to change two 
neutrons into two protons, it means that the one-body operator of H plays a less important role in the evaluation of 
M"q T , Second, in the case of the QBA the commutator [A(t/2), A(-t/2)] is just a constant but within the OEM it is 
a two-body operator. It means that the OEM approximations go beyond the QBA and RQBA |7^ , |75|] . 

Further, it has been shown that the energy difference between the initial and final nuclear states plays an important 
role in the der ivati on of the OEM transition operators and that the consideration of Coulomb interaction is crucial 
for th e OEM [ 1 74 1 . If the approximate nuclear two-body Hamiltonian H doe s not con tain a Coulomb interaction 
term [171,172 175 or the Coulomb interaction is considered in the way of refs. Irrcf 
transition operators is not consistent. Therefore 
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the derivation of the OEM 
an effective Coulomb interaction term Vc has to be considered, 



which represents partially the one body terms of the nuclear Hamiltonian 174 74.75]. The approximate OEM two- 
body nuclear Hamiltonian H containing Vc, central Vcn and tensor Vtn interactions is given as follows (the notation 
of ref. ]p77| is used): 



H « Vc + Vcn + Vtn, 



where 



Vc = lj2(E f - Ek) Q£ = j(l + r°)(l + Tf), 

V CN = [ ( SsE(n 3 ) KM) + 9so(r l3 ) K(ij) ) ^(ij) + 

( g TE {r l3 ) Il r e (ij) + g TO (r l3 ) U r a (ij) ) Hf (ij) }, 



(78) 



(79) 



V tn = ( 9tne { r ij ) K(*j) +9TNo(rij) K(ij) ) <%■ 



(80) 



(81) 



In the framework of the OEM2 for the nuclear matrix element M^j, one obtains 0,^5) 
M 2 GT = < 0;i^£ r ^ r / ( V sin3let (r i3 ) + 

V tri P let {rij) n?{lj) + ytensor^ g.. ) | Q + >5 



(82) 



where, 



—2 4 

ysmglet = ( g3 ) 

9te — gsE — ^qtne + A grE — gsE + ^9tne + A 
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9 so — 9to + ^9tno + A gso — 9to — %9tno + A ' ^ ^ 



11 1 2 

^tensor r ~ 



3 A -&9TNO + A &9tno + A 

1 1 



9so — 9to + ^9tno + A gso ~ 9to — ^9tno + A 



(85) 



Here V denotes the Principle value integration. We note that the central and tensor interactions appear only in the 
denominators of the OEM potentials. As A (half energy release for 2i//3/3-decay process) is about few MeV, only the 
long range part of nucleon-nucleon potential is relevant for the calculation. Therefore, the OEM potentials are not 
expected to be sensitive to the type of realistic effective interaction, e.g. Bonn or Paris ones. 

For the calculation of M"q T in Eqs. fl82] ) and (|8a) it is necessary to know the wave functions of the initial and final 
nucleus. The OEM c an be combin ed with the ground state wave functions of the QRPA or renormalized RQRPA 
models (OEM+RPA) [fL 75 177 17S ] . Recently, the consistently derived OEM transition operators have been c omb ined 



with pn- RQRPA wave functions of the initial and final nuclear systems to calculate 2z//3/3-decay of 76 Ge [17S]. A 
strong dependence of the matrix element Mq^ has been found on the choice of the pn-RQRPA ground state wave 
functions. Therefore, there is an interest to combine OEM-transition operators in Eqs. fls^ ) and ( p5|) with other 
ground state wave functions, e.g. shell model ones. 



D. Alternative methods 



Till now no consistent many body approach is available for the calculation of the many-body Green function 
governing the 2^/3/3-decay process, because of the computational complexity of the problem. The small predictive 
power of the QRPA is the motivation that stimulated us to seek an alternative way, which might produce more reliable 
results. 

The most promising are those methods, which can guarantee the particle number conservation and an explicit 
consideration of the Pauli principle. The two-vacua RPA (TVRPA) fulfill both criteria [180]. The trial nuclear wave 
functions of the intermediate odd-odd nucleus are constructed as linear combinations of proton-particle neutron-hole 
excitations from the ground state of the initial even-even (N,Z) nucleus and at the same time, as neutron-particle 
proton-hole excitations from the ground state of the final (N-2,Z+2) nucleus. Till now, the method has been applied 
only to the 2i//3/3-decay of 48 Ca. In order to understand its predictive power, it is clearly necessary to perform more 
calculations of the beta strength distributions and double beta decay matrix elements. 

In the past few years a new method for the solution of the nuclear many-body problem has been developed, 
which uses the shell model Monte Carlo technique (SMMC) [184]. The SMMC method enforces the Pauli principle 
exactly and scales more gently with the problem of the single particle basis size than the traditional shell model 
direct-diagonalization approach. The SMMC has been found successful in description of the Gamow- Teller strength 
distributions and has been also a pplie d for the 2i//3/3-decay of 48 Ca (complete pf shell) and 76 Ge (0f 5 / 2 ,lp,0g 9 /2 model 
space) giving promising results [184]. However, the dependence of M"q T on effective interaction and single particle 
energies remains to be investigated. 

Other approaches looking for the physical cause of the suppression mechanism of the 2^/3/3-decay are based on the 
approximate Wigner spin-isospin SU(A) aT symmetry in nuclei 1 181 182f| . The possibility of using SU(A) aT symmetry 



to analyze nuclear excitations (Gamow- Teller and isobaric analog resonances) has been discussed e.g. in [ 183 1 . We 
know that SU(A) aT is not an exact symmetry and is broken by the spin-orbit and Coulomb forces. It is worthwhile to 
notice that if the exact SU(4) aT symmetry is considered, 2i//3/3-decay is forbidden |8l|) , i.e. 2^/3/3-decay is determined 
by the mechanism responsible for breaking of this symmetry. In the approach proposed in ]182| the spin-orbit part of 
the mean field of the nucleus was assumed to be responsible for the symmetry breaking. 



V. RESULTS AND DISCUSSION 



A. QRPA and violation of the Pauli principle 

Most of the calculations of the 2z//3/3-decay matrix elements for nuclei heavier than 48 Ca rely on QRPA. However, 
the calculated matrix elements are uncertain because of their great sensitivity toJ=l + ,T = particle-particle 
component of the residual interaction governed by the g pp parameter ]63|-|65|]. The 2^/3/3-decay matrix element M"q T 
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is reduced by the enhanced ground state correlations with increasing strength of proton-neutron interaction and it 
even crosses zero. The dependence of M"q T on g pp for the 2vf3{3-dec&y of 76 Ge calculated within the pn-QRPA for 
three different model spaces is presented in the Fig. (||)a Jt5| . For 9, 12 and 21 level model spaces the collapse of 
the pn-QRPA occurs near the critical interaction strength given by g pp values 1.3, 1,1 and 0.9, respectively. A larger 
model space means more ground state correlations, i.e. a collapse of the pn-QRPA solution for smaller g pp . 

It was long believed that this curious behavior of the 2^/3/3-decay matrix elements in the QRPA is typical only for 
this decay mode. Recently, it has been found that in the framework of the pn-QRPA for a large enough model space 
the 0z//3/3-decay matrix elements shows similar instability with respect to g pp as the one known from the 2^/3/3-decay 
mode. The value of the matrix element crosses zero and it is then difficult to make definite predictions about the 
transition rate [see Fig. (Jjj)a]. This fact is not surprising as the origin of this dependence is the product of two 
one body transition densities of the charge changing operator to the lowest (or one of the lowest) 1 + state of the 
intermediate nucleus namely 

< o; || [c+X']i II It/ >< lt/|lti >< lt< II Kk^ih II 0+ >, (mi =m f = 1), (86) 

which is present in the calculation of each the double beta decay matrix element via the QRPA and its extensions (see 
Eqs. |f77f ) except the simple double Fermi transition. The theoretical nuclear matrix elements show large uncertainties, 
since they depend on the ground state correlations, which are overestimated in QRPA. 

It is worthwhile to notice that there is no mechanism in QRPA, which prevents the increase of the ground state 
correlation to unreasonable large values with increasing g pp , leading finally to a collapse of QRPA solution. The 
investigation of the M"q T within the pn-RQRPA j70| and the full-RQRPA [jn] has shown that the main source of 
the ground state instability in QRPA is the violation of the Pauli exclusion principle, if one has large ground state 
correlations. 

The RQRPA, which considers the Pauli exclusion principle in an approximate way, shifts the collapse of the QRPA 



outside the physical range of g pp and shows a less sensitive dependence of M"q T on g pp |70]-|72|,185|. Fig. ([|)b shows 
Mq 1 ^ for the 2z//3/3-decay of 76 Ge calculated by pn-RQRPA |75|. We see that the crossing point through zero is shifted 



to larger g pp in comparison with Fig. (|])a. A steeper slope of M"q T as a function of g pp for larger model spaces may 
indicate that the approximate treatment of the Pauli exclusion principle by the pn-RQRPA is for larger g pp not good 
enough. 

The RQRPA offers considerable less sensitive results to g pp in comparison with the QRPA also for the 2^/3/3-decay 
transitions to the + and 2 + excited states of the final nucleus Jl^,|l§]. The same is valid also for the calculatio n of 



M% aaa = M^ T (l - X f) associated with the neutr ino mass mechanism of the Oi//3/3-decay [see Fig. (g)] |109| , |l70| 

In both QRPA and RQRPA treatments of the double beta decay matrix elements one has two sets for the same 
intermediate nuclear states generated from initial and final nuclei, which do not coincide with each other. The overlap 
factor of these two sets is introduced artificially and the effect of this inaccuracy is not well understood. One can learn 
more about this problem, if the mismatching between the initial and final QRPA (RQRPA) Hamiltonians is studied 
J75|. It has been found that pn-RQRPA Hamiltonians demonstrate a better mutual agreement than the pn-QRPA 
ones |fr5fl . 

Recently, the validity of the RQRPA have been discussed within schematic models, whic h are try ing to simulate 
the realistic cases either by analytical solutions or by a minimal computational effort |75, 186- 188)] . Again it was 



confirmed that the renormalized QRPA is superior to the QRPA. But differences were found between the exact and 
the RQRPA solutions after the point of collapse of QRPA. However, it is not clear whether these results would also 
hold for realistic calculations with large model spaces, realistic effective NN- interactions and by consideration of the 
two- vacua problem. In fact we do not know of any exactly solvable realistic model. 

We note that there is a principal difference between the 2z//3/3-decay and Of /3/3-decay modes from the nuclear physics 
point of view. In the 2^/3/3-decay the correlation of two single nucleon beta decays in the nucleus are taken into account 
by the interaction part of the nuclear Hamiltonian. For a correct treatment of these correlations higher order terms 
in the boson (renormalized boson) expansion of the nuclear Hamiltonian are important [J75[. Several procedures 



have been proposed which outline the importance of the higher order effects J67,189,190, 188| . It is expected that 
the consideration of these higher order approximations will improve the correspondence between the initial and final 
nuclear Hamiltonian and more reliable results can be obtained. 

The Oi//3/3-decay is not a higher order effect in the boson expansion of the nuclear Hamiltonian as the dominant 
correlation of the two beta decay vertices in the nucleus comes from the exchange of lepton number violating par- 
ticles. Therefore the renormalized QBA scheme should be sufficient. We believe that the renormalized- QRPA is at 
present the most reliable method to deduce the desired interesting lepton number non-conserving parameters from 
the experimental lower limits on the half-lifcs of 0^/3/3-decay of heavier nuclei. 
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B. Limits on lepton number violating parameters 



It is a common practice to impose limits on different lepton number non-conserving parameters to assume that 
only one mechanism dominates at a time. A more accurate procedure requires the construction of multi-dimensional 
plots. This is worth doing only when O^/3/3-decay has been found in the future. 

In order to extract the limits on the effective neutrino mass < m„ > and effective Majoron coupling < g > from 
the experimental limits of O^/3/3-decay lifes the nuclear matrix elements M*q T and \f have to be calculated. We have 
evaluated these matrix elements for the experimentally most interesting A=76, 82, 96, 100, 116, 128, 130, 136 and 150 
nuclear systems within the framework of pn-RQRPA [|70|j7l|] . The same model space and nuclear structure parameters 



have been used as those considered in [153]. The calculated nuclear matrix elements for g pp — 1.0 together with upper 
bounds on < m u > and < g > are presented in Table IV. We see that nowadays 76 Ge provides the most stringent 



limit on the effect ive n eutrino mass: < m v >< 0.54eU. This value is about twice smaller as the value obtained within 
the full-RQRPA (|l09| , which include proton-neutron pairing correlations. The most stringent upper constraint on 
< g > is associated with A=128 and is equal to 5.5 x 10~ 5 . 

As w e mentio ne d alr eady many differ ent nuclear models (sh ell m odel ||l56| Jlj|, general seniority scheme [191], pn- 



QRPA |!67Ul68| ,fr9T|l26| ; full-QRPA |l2f| and pn-, full-RQRPA JTo| and others ) have been used for the calculation 
of the 0^/3/3-decay matrix elements. Unfortunately, there is no possibility to assign an uncertainty to the calculated 
theoretical values. One can only claim that the renormalized QRPA predicts more reliable results in respect to the 
QRPA because of reasons discussed before. A comparison between the shell model and the QRPA (RQRPA) is not 
possible as both models work with fully different model spaces and nucleon-nucleon interactions. We remind that the 
QRPA is loosing sense for a small model space and shell model studies for heavier open-shell nuclei are not possible. 
Thus one can expect that for the A=48 system the shell model results are more reliable and for other nuclear systems 
the RQRPA perhaps have more predictive power. 

First, we discuss here the light Majorana neutrino exchange mechanisms. The 0^/3/3-decay half-lifes for several 
nuclei calculated from some representative nuclear matrix elements of different nuclear models are listed in Table 
[v|. We have assumed that only one mechanism at a time is important and that < m v >= leV, < A >= 10~ 6 , 
< rj >= 10~ 8 and < g >= 10~ 5 . By glancing at Table [v| we see that most theoretical works have concentrated on 
A=76, 82, 128 and 130 nuclear systems and that the calculated half-lifes agree within a factor ten. It is remarkable 
that there is a large disagreement between the two different shell model calculations for 76 Ge and 82 Se Ol56^ . Wo 
note that from the theoretical point of view the best candidates for the experimental measurements are ^ u 7Vd and 
130 Te nuclei with the shortest half-lifes. The upper bound on < m v >, < A >, < r\ > and < g > can be deduced 



from the experimental limit on the 0i//3/3-decay given in Table | with help of the theoretical half- life T*^ or [ see 



Table as follows: 



</> = Ci 



\ 



rptheor 
± X/2 

rpCXp • 

1 l/2 



(87) 



where Gj takes the value 1 eV, 10 6 , 
limits are as follows: 



10 8 and 10 for I = m„, A, r\ and g, respectively. Thus, the most restrictive 



<m v > < 0.4-1.3eV [ 7fa Ge, ref. 21], 



0.7- 1.5eV [ 128 Te, ref. 35] 



< A > < (0.8 - 2.1) x 10~ 6 [ 76 Ge, ref. 21] (2.5 - 7.3) x 10~ 6 [ 136 Ae, ref. 37] 
<?/> < (0.4- 1.8) x 10" 8 [ 76 Ge, ref. 21], (0.9 - 3.7) x 10" 8 [ 128 Te, ref. 35] 

< g > < (2.9 - 5.8) x 10~ 5 [ 128 Te, ref. 35], (9.8 - 22.) x 10~ 5 [ 100 Aio, ref. 29] 



(89) 
(90) 
(91) 



The sensitivity of different experiments to these < m„ >, < A >, < r\ > and < g >. is drawn in Fig. ^. Currently, 
the Heidelberg-Moscow experiment pfl| offers the most stringent limit for < m v >, < A >, < 77 > and the 128 Te 
experiment |35| for < g >. 

Second, we present nuclear matrix elements associated with the Ifip SUSY and heavy neutrino exchange mechanisms. 

The nuclear matrix elements governing the two-nucleon, one pion-exchange and two-pion exchange J]L V SUSY 
contributions to the Oi^/3/3-decay have been calculated only within the pn-QRPA, pn-RQRPA and the full-RQRPA 
approaches |58|,|59,153]. It has been found that the pion-exchange contribution dominates over the conventional two- 
nucleon one and its contribution can be safely neglected J59[ . The calculated pion-exchange nuclear matrix elements 
^ttTV f Qr £ ne 0i//3/3-decay of various isotopes within the pn-RQRPA (g pp = 1.0) are presented in Table VI [ 153 1 . 



The two largest are -1073 and -756 for A=150 and 100, respectively. It has been found that the variation of M.'' 



-v 
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on g pp do not exceed 15% within the physical region of g pp . In the evaluation of M* N the pion vertex structure 
coefficient was considered to be a 27r — 0.2 (VIA) ]59[| , We note that the value of M* N calculated for 76 Ge within the 
full-RQRPA [5S| ] is suppressed by about 30% in comparison with pn-RQRPA. The upper bounds on r)^ av [see Eq. 

we see 



(p9|)] have been deduced by using the most stringent lower limits of T^J£ m Table |. By glancing at Table VI 



that 76 Ge provides the most stringent limit: r]f£P — 5.5 x 10 9 . By using the Eqs. © and © one obtains for the 
1st generation Ifi p Yukawa coupling constant |l53| 



A' ni < 9.3 x 10- V m * Y( m * ) 1/2 . (93) 
111 - VlOOGeV/ VIOOGeV/ v ' 



These limits are much stronger than those previously known and lie beyond the reach of the near future accelerator 
experiments. The size of A' ln depends on the masses of the SUSY-partners. If the values of these masses would be 
around their present experimental lower limits ~ lOOGeV [Q, one could constrain the coupling to A' 1X1 < 1.3 ■ 10~ 4 . 
A conservative bound can be set by assuming all the SUSY-masses being at the " SUSY- naturalness" bound of 1 TeV, 
leading to A' m < 0.04. 

The nuclear matrix element M. n associated with the exchange of heavy Majorana neutrino are shown for different 



nuclear systems in Table VI . They have been evaluated within the pn-RQRPA for the same nuclear model parameters 
as Mq^, and \f hi Table IV. The two largest elements are -206 and -151 for for A=150 and 100, respectively. It 



is worthwhile to note that these values can be reduced, if higher order terms of the nucleon current are considered 



[ 192 1 . For the A=76 system the most stringent limit on the relevant lepton number non-conserving parameter has 
been found: 77^ < 2.6 x 10 -8 . 

VI. CONCLUSION AND OUTLOOK 

In this review we have discussed the recent progress in the field of the nuclear double beta decay and its present 
status. 

The 2^/3/3-decay remains at the forefront of nuclear physics. The established 2^/3/3-decay half-lifes for a couple of 
isotopes constrain nuclear theory and stimulate its further development. There are still a number of open questions 
concerning the solution of the nuclear many-body problem. In order to increase the predictive power of the nuclear 
many-body models it seems to be necessary to include better the Pauli exclusion principle in the evaluation of the 
2j//3/3-decay matrix elements. This is done in the renormalized QRPA. It is an improvement over the commonly used 
QRPA approach. An improved description of the complete set of the intermediate nuclear states, which are needed 
for the calculation of double beta decay, can possibly also be obtained with growing computer power for example in 
the shell model approach. 

Additional experimental information about the 2j//3/3-decay and related processes is of great interest. We note 
that the NEMO experiment allows to perform a precise measurement also of the energy and angular distributions 
of the outgoing electrons. Further, there is a chance to detect the 2j//3/?-decay to the excited + and 2+ states of 
the final nucleus at the level of 10 21 -10 22 years with present low-background detectors. New (3(5 experiments are 
planned by the detection of the capture of two bound atomic electrons ]193[ . A recent theoretical analysis has shown 
that the experimental measurement of the double beta decay induced in a neutrino beam can be feasible due to an 



enhancement due to a resonance in the corre sponding amplitude [194]. Such an experiment is planned for the double 



beta decay at the Kurchatow Institute | 194 |. These future experiments are important because they complement the 
known informations from the conventional 2i//3/3-decay. 

The 0^/3/3-decay is a powerful tool for the study of lepton number conservation in general. In particular, it 
contributes to the search for the Majorana neutrino mass, the heavy W-boson (mainly responsible for a possible 
right-handed weak interaction), the J^ p SUSY and leptoquarks on the sensitivity level competitive to Large Hadron 
Colliders. We remind that the study of rare decays at low energy is complementary to high energy physics experiments. 

Searches for double beta decay are pursued actively for different nuclear isotopes, but the O^/3/3-decay has been 
not yet observed, experimentally. The currently laboratory upper limits on the half-life provide severe constraints on 
the effective Majorana neutrino mass < m„ >, on parameters of the left-right symmetric models < A > and < 77 >, 
on the effective Majoron coupling constant < g >, on the i?-parity violating coupling constant A' ln and on the the 
inverse effective mass of heavy neutrinos 77^ as follows: 
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<m v > < 0.4-1.3eV, <g> < (2.9 - 5.8) x 

<A> < (0.8-2.1) x <r;> < (0.4 - 1.8) x 

V >" ^ 1 ^ m A^)\^hj) 1 ''- •»< 2.6X10- (94, 

The uncertainty of these parameters is due to the ambiguity of nuclear matrix elements. 

The above limits constitute our most stringent test of the lepton number conservation. It is worthwhile to notice 
that the Oi//3/3-decay imposes very restrictive bounds on the lepton number violating sector of R-parity violating 
SUSY models. The upper limit on A' m (deduced from the pion-exchange SUSY mode for the O^/3/3-decay |3j|), is 
significantly stronger than those previously known or expected from the ongoing accelerator experiments. The §v(3(3- 
decay constraints on lepton number violating parameters in Eq. (^) must be taken into account by theoreticians, 
when they build new theories of Grand Unification. 

In the 0^/3/Ldecay it is still possible to improve the existing limits. The extension of the ongoing experiments: ELE- 



GANT ( 48 Ca) |195| , Heidelberg-Moscow (™Ge) gj, NEMO ( 100 Alo, 116 Cd) f^Q, the 130 Te cryogenic experiment 
p6| and Caltech Neuchatel TPC 136 Ye |37) is expected to reach the level of 0.1 eV. The p rospect for exploration of 
SUSY in the next generation of the 0!//3/Ldecay experiments have been discussed in [ |196| . Recently, a new project 
GENIUS for the measurement of the Oi^/3/3-decay of 76 Ge have been proposed, which suppose to probe neutrino masses 



down to 10 2 -10 3 eV by using 1 ton of enriched 76 Ge [ 197 1 . If the O^/3/3-decay will be discovered, it would be a major 



achievement, which would initiate new experimental activities both for rare decays at low-energies and also at high 
energy accelerator facilities. In the case that the lepton number violation is out of reach of near future Oz^/3/3-decay 
experiments, the improved upper limits will yield more stringent constraints for the Grand Unified Theories and Super 
Symmetric (SUSY) models. 
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TABLE I. Experimental results for 2vj3f3-, Ov/3/3- and 0z//?/3</>-decay modes for A=48, 76, 82, 96, 100, 116, 128, 130, 136, 150, 
238 and 242 nuclear systems. gch.(rch.) - geochemical (radiochemical) data. 
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TABLE II. The transition operators 0j(rj, r,-, <Ti, Oj) = const- ^ Vi(q,r)dq <S/(f^, r+ij, (Ji, <7j) (see Eqs. (p[)) for 
the nuclear matrix elements associated with exchange of light and heavy particles are listed ( = r,: - r-j, r+ij — rt + r-j, 
fy = r ij'/l r ij|> r+ij = r+ij/|r+ij|, r = |fijj ). The following notations are used: R is nuclear radius, ji (i=l,2) are spherical 
Bessel functions. f(q 2 ) = 1/(1 + q 2 /m 2 A ) 2 is the dipole nucleon form factor with cut-off mA=850 MeV. and m p are the 
pion and nucleon masses, respectively. A is the closure energy, which can be found for different nuclei e.g. in [79]. Sij =3o~i ■ 
iij (7 j ■ iij - Gi ■ (7 j is the tensor operator. <r_ij = Oi — Oj. gy = 1.0, qa = 1-24 and [ip = 4.7 
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TABLE III. The nucleon structure coefficients of the two-nucleon mode for the R v SUSY mechanism for the O^/3/3-decay 
calculated within the non-relativistic quark (QM) and Bag models. They have been calculated using Table I and Eqs. (54-57) 
ofref. [58]. 
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TABLE IV. The matrix elements Mg T and X f of the Oi//3/3-decay for 48 C*a, 76 Ge, 82 Se, 9(i Zr, 100 Mo, 116 Cd, 128 Te, 130 Te and 
136 Xe calculated in the framework of pn-RQRPA. The upper limits on the lepton number non-conserving parameter < m„ > 
(< g >) have been extracted from the best presently available experimental Of/3/3-decay (O^/3/30-decay) half-lifes [see Table I ]. 
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TABLE V. Calculated O^/3/3-decay half-lifes for < m v >= leV, < X >= 10~ 6 , < r? >= 10" 8 and < g >= 10" 5 
that one mechanism dominates at a time. The results have been obtained within the following nuclear models: 
[1,156], generalized seniority scheme [191], quark confinement model + pn-QRPA [132], pn-QRPA [79,167,168] 
full-QRPA [126]b, pn-RQRPA (p.w. - present work) and full-RQRPA [109]. 
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TABLE VI. The nuclear matrix elements for SUSY (.M ) and heavy neutrino exchange (Mn) 0^/3/3-decay (see Table 
0) for A=76, 82,96, 100, 116, 128, 130, 136 and 150 nuclear systems calculated within pn-RQRPA [153]. The lepton number 
non-conserving parameters T] SU3Y and t]n have been deduced from the most stringent experimental limits of 0^/3/3-decay lifes 
[see Table [jj]. The upper bounds on SUSY coupling constant Am have been obtained by using rriq = m g = = 100 GeV [ see 
Eqs. ( |60[ ) and (|6l[)]. T®J~ suay anc [ rpOj-heavy arg ^ e calculated half lifes assuming r/ SUSY — 10~ 9 and t)n = 10~ 8 , respectively. 
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FIG. 1. Mechanisms for the O^/3/3-decay associated with the exchange of light Majorana neutrino: (a) the light neutrino 
mass mechanism, (b) the Majoron mechanism, (c-d) right-handed current mechanisms. The following notation is used: Wl(h), 
(1l(r) and &l(r) are left- (right-) handed u-quark, d-quark and electron, respectively. W is vector boson (light or hypothetical 
heavy) and i>i (i=l,2...) is the Majorana neutrino). 

FIG. 2. Feynman graphs for the supersymmetric contributions to O^/3/3-decay. Ul, cLr and ei, have the same meaning as 
in Fig. [[]. ul, d.R and El are left-handed u-squark, right-handed d-squark and left-handed selectron, respectively. % and g are 
neutralinos and gauginos, respectively. 

FIG. 3. Two-nucleon (a), one-pion (b) and two-pion (c) mode of neutrinoless double beta decay in imparity violating SUSY 
extensions of the standard model. 



FIG. 4. The nuclear matrix element Mgt for the 2i//3/3-decay of 76 Ge calculated within the pn-QRPA (a) and pn-RQRPA 
(b) is plotted as a function of the particle-particle coupling constant g pp . The dashed line corresponds to the 9-level model 
space (the full 3 — 4hui major oscillator shells), the dot-dashed line to the 12-level model space (the full 2 — 4hu) major oscillator 
shells) and the solid line to the 21-level model space (the full — 5hui major oscillator shells). 



FIG. 5. The calculated nuclear matrix element M% aaa = Mq t (1 - xf) for the 0^/3-decay of 76 Ge, 100 Mo, 116 Cd, 128 Te 
and 136 Xe as a function of the particle-particle interaction strength g pp . In (a) M^ aaa has been calculated with the pn-QRPA 
method, in (b) with the full-RQRPA . 

FIG. 6. The sensitivity of different experiments to the lepton- number violating parameters < m v >, <A>, < r\ >, < g >, 
Vsusy an d VN- In construction of the histograms the best available limits on Qvj3j3- and O^/3/30-decay half-lifes from Table I 
have been used. The open and black bars in (a-d) correspond to the smallest and largest value of Tj/| or for a given nucleus in 
Table V, respectively. 
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